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ABSTRACT 


Employing  an  analogy  between  thermally  Induced  and  irradiation  induced 
creep,  physical  arguments  are  used  first  to  deduce  a one-dimensional 
constitutive  relation  for  metals  under  stress  in  a high  temperature  and 
high  neutron  flux  field.  This  constitutive  relation  contains  modified 
superposition  Integrals  in  which  the  temperature  and  flux  dependence  of 
the  material  parameters  is  included  via  the  use  of  two  reduced  time  scales; 
linear  elastic,  thermal  expansion  and  swelling  terms  are  also  Included. 

A systematic  development  based  on  thermodynamics,  with  the  stress,  tempera- 
ture increment  and  defect  density  increment  as  independent  variables  in  the 
Gibbs  free  energy,  is  then  employed  to  obtain  general  three-dimensional 
memory  integrals  for  strain;  the  entropy  and  coupled  energy  equation  are 
also  obtained.  Modified  superposition  Integrals  similar  to  those  previously 
obtained  by  physical  argument  are  then  obtained  by  substituting  special 
functions  into  the  results  of  the  thermodynamic  analysis,  and  the  special 
case  of  an  isotropic  stress  power  law  is  examined  in  detail. 


1 . INTRODUCTION 


Metals  under  stress  at  elevated  temperatures  experience  time  dependent 
strain,  i.e.,  the  well-known  phenomenon  called  thermally  induced  creep.  As 
general  rule,  thermally  induced  creep  in  a metal  is  considered  Insignificant 
at  temperatures  below  one  third  of  the  melting  point.  Similarly,  metals 

under  stress  in  a high  neutron  flux  field,  even  at  temperatures  well  below 
one  third  of  the  melting  points,  also  experience  significant  creep  strains. 

This  less  well-known  phenomenon  is  called  irradiation  induced  creep.  Gilbert 
[1]  summarized  some  of  the  possible  microscopic  mechanisms  associated  with 
irradiation  in  metals,  and  observed  that  irradiation  may  "enhance"  thermal 
creep  mechanisms  as  well  as  inducing  new  creep  mechanisms.  It  is  believed  that 
the  most  significant  irradiation  effect  in  a metal  is  the  production  of  point 
defects  by  the  collision  of  high  energy  neutrons  with  the  lattice  atoms. 

Thus  the  production  of  point  defects  forms  an  essential  element  in  any  study 
of  the  enhancement  of  thermal  creep  and/or  the  inducement  of  irradiation  creep. 

The  main  purpose  of  this  paper  is  to  construct  material  constitutive 
relations  valid  in  a neutron  flux  field  at  high  temperatures.  Since  Andrade's 
initial  work  in  thermal  creep  early  in  this  century,  many  mathematical  models 
have  been  proposed  to  represent  a great  accumulation  of  constant  stress  data. 
These  models  can  be  generalized  to  variable  stress  by  using  the  strain  or 
time  hardening  techniques  (see  [2]),  and  reasonably  good  results  may  be 

obtained.  In  [3],  a single  integral  form  of  the  constitutive  relation  is 
obtained  by  applying  the  strain  hardening  technique  to  each  thermal  creep 
component,  and  we  follow  this  approach  in  the  present  paper.  Also,  since  the 
thermal  creep  properties  depend  strongly  on  temperature,  a thermorheologically 
simple  material  Is  assumed  and  a thermal  reduced  time  is  employed  [4]  in  the 


thermal  creep  compliance  function. 

As  to  irradiation  induced  creep,  a limited  amount  of  data  has  been 
obtained  during  the  past  two  decades  [5  - 9]  for  various  metals  and  alloys 
in  various  geometries  (e.g.  rods,  spring,  etc.).  Such  data  are  generally 
for  the  case  of  constant  flux  and  constant  stress,  and  it  is  found  that  the 
creep  rates  are  dependent  on  the  materials,  their  pre-irradiation  treat- 
ments, the  temperature  level,  the  stress  level  and  the  flux  level.  A 
reasonable  fit  for  these  data  is  obtained  by  a creep  model  characterized 
by  the  sum  of  linear  elastic,  steady  creep  and  transient  creep  strain  com- 
ponents. The  strain  hardening  technique  may  be  applied  to  generalize  this 
type  model  to  the  case  of  variable  stress, which  may  then  be  expressed  in 
terms  of  memory  integrals.  The  flux  dependence  of  the  material  parameters 
in  irradiation  induced  creep  plays  a role  similar  to  the  temperature 
dependence  of  the  parameters  in  thermally  induced  creep,  and  thus  we  intro- 
duce the  concept  of  a "fluxorheologically^simple  material  and  use  an 
irradiation  reduced  time  in  the  irradiation  creep  compliance  function.  At 
high  levels  of  both  temperature  T and  flux  <f>,  thermal  and  irradiation 
creep  occur  simultaneously  (i.e.  "thermoirradiation"  induced  creep),  with 
the  material  parameters  depending  on  both  T and  <f>.  For  this  case,  Rashid 
[10]  proposed  a single  memory  integral  with  a modified  reduced  time  contain- 
ing functions  of  T and  $.  In  this  paper  we  use  an  alternate  approach 
whereby  we  superpose  memory  integrals  with  two  distinct  reduced  times,  where 
both  times  contain  functions  of  T and  $ with  one  dominated  by  the  thermal 
dependence  and  the  other  by  the  flux  dependence. 

Constitutive  relations  are  state  equations  in  thermodynamics,  and  thus 
they  must  be  consistent  wlth^first  and  second  laws.  Thermodynamic  theory 
has  been  applied  to  viscoelastic  materials  along  basically  two  different 
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approaches:  (1)  The  energy  functional  approach  of  Coleman  [11],  and  (2) 

The  state  variable  approach  of  Biot  [12]  and  Schapery  [13].  Similar 
results  may  be  obtained  -ffcrough  either  approach,  as  discussed  in  [14]. 

Chang  and  Cozzarelli  [15]  employed  the  energy  functional  approach  to  obtain 
a single  integral  constitutive  relation  for  nonlinear  thermoviscoelastic 
materials.  A similar  approach  was  used  in  [16]  to  develop  constitutive 
relations  for  thermally  and  irradiation  induced  creep;  the  present  paper  is 
essentially  an  amplification  and  generalization  of  that  initial  effort. 

Although  the  neutron  flux  $ in  irradiation  creep  produces  an  effect 
analogous  to  that  of  the  temperature  T in  thermal  creep,  $ is  an  external 
input  and  not  an  internal  thermodynamic  variable  as  is  T.  However,  Makin, 
Whapham  and  Minter  [17]  found  that  if  annealing  is  neglected  the  rate  of  change 
of  the  imperfection  density  n is  proportional  to  the  flux  <f>.  Based  on 
this  observation,  we  employ  the  imperfection  density  n instead  of  the 
neutron  flux  <t>  in  our  thermodynamic  development. 

An  important  contribution  to  the  total  strain  in  a high  radiation  field 
is  due  to  the  swelling,  which  is  analogous  to  the  thermal  expansion  in  a 
thermal  field.  Swelling  is  defined  as  the  volume  change  due  to  neutron  flux 
in  the  absence  of  applied  stress.  This  phenomenon  was  discussed  some  time 
ago  by  Cawthorne  and  Fulton  [18],  and  since  then  various  empericai  equations 
have  been  developed.  In  this  paper,  we  assume  that  the  swelling  is  simply 
proportional  to  the  fluence  $t,  as  in  [19,  20].  Since  d>  has  been  taken 
here  proportional  to  A, the  swelling  is  thus  proportional  to  the  increment 
in  defect  density.  This  establishes  a direct  analogy  between  swelling  and 
thermal  expansion,  which  is  of  course  proportional  to  the  Increment  in 
temperature. 

In  Section  2,  simple  creep  models  deduced  from  one-dimensional  creep 
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test  data  are  generalized  to  the  case  of  variable  stress,  temperature  and 
flux,  resulting  in  modified  superposition  integrals.  In  Section  3,  the 
thermodynamic  development  of  the  general  three-dimensional  single  integral 
form  of  the  constitutive  relation  is  discussed  in  detail.  In  Section  4 , 
we  show  that  the  superposition  integral  representation  in  Section  2 is  in 
fact  a special  case  of  the  general  representation  in  Section  3,  obtainable 
by  selecting  a special  form  for  the  Gibbs  free  energy.  In  Section  5,  one- 
dimensional expressions  and  compressibility  effects  are  examined  for  the 
special  case  of  isotropic  stress  power  laws.  A short  summary  then  follows. 
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2.  PHYSICAL  PRELIMINARIES 

In-reactor  creep  strain  c in  metals  is  a complex  function  of  stress 
o,  temperature  T,  and  flux  <t>.  For  convenience,  we  classify  it  into  three 
categories:  (1)  Thermally  induced  (high  T,  low  4>> , (2)  Irradiation 

induced  (low  T,  high  d>> , (3)  Thermoirradiation  induced  (high  T,  high  <f>). 

One-dimensional  constitutive  relations  for  each  of  these  categories  of  creep 
are  developed  here  by  various  physical  arguments  based  on  experimental 
observation. 

2(a).  Thermally  induced  creep 

Many  models  for  thermal  creep  are  given  in  [2],  resulting  in  a variety 

of  strain-time  relations  which  contain  the  effects  of  temperature  and  stress. 

Following  the  approach  in  [3],  we  write  the  strain  e,  at  time  t,  for 

constant  temperature  T and  constant  stress  r , as  a sum  of  linear 

o o 

elastic,  thermally  induced  steady  creep,  thermally  induced  transient  creep, 
and  thermal  expansion  components,  i.e., 

t « + £'’<«;>  #r‘T.>+  F^‘(<r.)  [ i - * otTeo  (i) 

Here,  E is  the  elastic  modulus,  F are  functions  of  stress,  ^ is 

a function  of  temperature,  A ^ are  constants,  and  is  the  coefficient 

of  thermal  expansion  for  temperature  increment,  0q  ■ T0'*T’ri  relative  4r» 
constant  reference  temperature  T^.  In  the  remainder  of  this  section,  we 
first  generalize  Equation  (1)  to  the  case  of  variable  stress  and  then  generalize 
it  to  the  case  of  variable  temperature. 

Generalization  to  Variable  Stress  d’(t) 

The  first  term  in  Equation  (1),  the  elastic  strain,  is,  of  course, 
already  valid  for  variable  stress  d(t).  Various  techniques  may  be  employed 
to  extend  the  second  and  third  terms  to  the  case  of  variable  tf(t).  Follow- 
ing the  strain  hardening  approach  in  [3],  we  assume  that  each  component  of 
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.00 


creep  strain  c may  be  written  as  a function  of  stress  and  strain,  i.e., 

(2) 


<*) 

|f  = i [<r(t),£tW<4» 


at 

Applying  this  approach  to  Equation  (1)  term  by  term,  we  obtain 


ai«<t)  ci) 


at 


(3) 


i = 2,--,M  (4) 

where  tg  is  the  steady  creep  and  e is  the  i-th  component  of 

transient  creep.  With  the  use  of  an  integration  factor  and  an  integration 
by  parts,  an  integral  form  of  Equation  (1)  valid  for  variable  stress  may  be 
obtained  as 


t = 


+ f #*)£(*>&'  * 1 a‘;V.)  exf>C-/£%(T0rt] 

o T 


• J expC  A,-  $r(.T,)t']  Ft  (<r)dt  ♦ aT0o 

O 

Equation  (5)  may  also  be  rewritten  in  terms  of  memory  integrals  as 


JO. 


(5) 


aF‘U(«r) 


t - -r  * J 5T(T.)(t-t')  — 


at' 


♦£  j*  | I - CXp  ‘t't 


3FtW)  . , 


+ 0<T  08  (6) 

The  integrals  in  Eq.  (6)  are  modified  superposition  Integrals  similar  to  the 
type  first  Introduced  by  Leaderman. 

Generalization  to  Variable  Temperature  T(t) 

We  assume  that  the  temperature  dependence  of  the  material  properties  is 
of  the  thermorheologlcally  simple  type.  As  in  [4],  we  Introduce  a reduced 
time  as 

£t  * £ #T<T(t')> 


(7) 
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and  transform  real  time  space  (x^  t)  Into  (x£,  CT>  space.  For  clarity,  we 
place  a carat  above  each  dependent  variable  In  the  new  space,  for  example 


£<Xi.$T>  5 f (Xi.t)  * -f  <X;,  3 


(8) 


where  g-A(x  ,£T)  is  the  inverse  of  Equation  (7  ). 
Accordingly,  Equation  (6)  becomes 

A t 'Ml)  a 

A <T  rT  * * * N B Ft  ( <T ) . / 

£ -X  * L 5r 

A (c)  /V 

M 4-  . <i>  , , , d Ft  (<r>  J 

*z, r «»p£-Ar <4r-5;»)  — «r" 

,«2  o ~T 


+ otT  9 

Finally,  Equation  (9)  may  be  rewritten  as 

s-4  ♦.£, 

where  JT^(5T>  are  the  thermal  steady  and  transient  creep  compliance 
functions  ( $T)  = and  JT  <^T)  = I - exp  (-At  £-r^»  ‘ = 


(9) 


(10) 


Equation  (10)  may  be  represented  by  a non-linear  generalized  Kelvin  model, 

and  the  terms  in  the  summation  are  modified  superposition  integrals  for  a 

thermorheologlcally  simple  material.  In  the  case  where  aging  effects  are 

significant,  a time  hardening  procedure  may  be  used  to  generate  creep  com- 
(i)  • 

pllances  JT  (^,  £^)  as  in  [3];  aging  effects  are  also  discussed  in 

detail  by  Krempl  [21]. 

2(b).  Irradiation  Induced  Creep 

For  emphasis,  we  summarize  here  some  pertinent  earlier  remarks  on  irradia- 


tion creep  given  in  the  introduction.  Some  data  has  been  obtained  for  steels 
and  alloys  by  various  constant  stress  and  flux  tests,  such  as  the  stretched 
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spring  [5],  pressurized  tube  [6],  bent  beam  [7],  torsion  rod  [8]  and  uniaxial 

tension  [9],  as  summarized  in  [19].  Although  reliable  data  are  still  scarce 

(especially  at  high  f luences) , a reasonable  model  can  be  constructed.  It  is 

found  that  the  expression  for  irradiation  induced  strain  at  constant  flux 

(<J>q)  and  constant  stress  (oq)  is  similar  to  Equation  (1).  However,  the 

material  creep  parameters  now  depend  on  a function  of  the  flux  rather  than  of 

the  temperature,  and  the  thermal  expansion  term  in  Equation  (1)  is  replaced 

by  a swelling  term.  It  is  also  found  that,  in  the  absence  of  aging  (i.e. 

annealing),  the  flux  is  proportional  to  the  rate  of  change  of  the  point 

defect  density  n.  For  the  convenience  of  our  thermodynamic  development, 

we  choose  to  express  the  material  creep  parameters  as  functions  of  n 

rather  that  $.  The  swelling  term  is  often  less  significant  than  the  other 

strain  terms,  and  we  express  it  simply  as  linear  in  <p  t.  Since  4>  is 

o o 

proportional  to  n^,  we  may  rewrite  this  term  as  linear  in  the  imperfection 
density  increment  * nQ-nr  relative  to  the  reference  n^.  Accordingly, 

£ = — + F aW  *„<  ry  t 


N l i)  , ti)  • 

fk  ( cr0)  C i - •+  olRyo 


(ID 


where  F ^ are  functions  of  stress,  £ is  a function  of  n (i.e.  flux), 

K K O 

A^^  are  constants  and  aR  is  the  coefficient  of  swelling.  Equation  (11) 
is  remarkably  similar  to  Equation  (1),  but  in  addition  to  the  replacement  of 


T by  n and  6 by  v , it  differs  in  two  important  respects:  (1) 
o o o o 

Whereas  the  F^^(oq)  are  highly  non-linear,  the  Fr^(oq)  are  often  chosen 
linear,  (2)  While  4>T  is  an  exponential-like  function  of  temperature,  jJ»R  is 
often  linear  In  n (i.e.  the  flux). 


often  linear  in  n 
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Proceeding  by  analogy  with  thermally  Induced  creep,  we  may  generalize 
each  irradiation  induced  creep  term  to  variable  stress  by  the  strain  hardening 
technique  (no  aging),  and  then  to  variable  n by  introducing  an  irradiation 
reduced  time 

i.-  (12) 

Thus , 


£ = 


aF„  (<r) 


+ 


A 

,v 


(13) 


where  are  the  irradiation  induced  creep  compliance  functions.  Continuing 

the  analogy  with  thermally  induced  creep,  we  designate  the  material  as 

fluxorheologically  simple.  Note  that,  for  all  F ^ (a)  linear,  the  summation 

K 

in  Equation  (13)  may  be  dropped  if  a summed  compliance  is  introduced.  Finally, 
we  note  that  aging  may  be  more  significant  in  irradiation  creep  than  in 
thermally  induced  creep,  and  thus  it  may  be  necessary  to  employ  a time  harden- 
ing procedure  to  generate  the  irradiation  creep  compliance  functions 

JR  (CR’  V’ 

2 (c).  Thermoirradiation  Induced  Creep 

At  high  temperature  and  high  neutron  flux,  the  constitutive  relation  will 
likely  be  some  complicated  extension  and  combination  of  Equations  (10)  and  (13). 
Very  little  is  known  about  this  relation  at  present,  but  it  appears  that  the 
thermally  induced  creep  parameters  will  now  also  depend  on  some  weak  function 
of  flux  jtT(n)  i.e.,  irradiation  enhancement).  For  example,  in  [22],  the 
thermally  induced  creep  rate  is  increased  a factor  equal  to  (1  + Bn). 
Similarly,  the  irradiation  induced  creep  parameters  will  depend  on  some  weak 
function  of  temperature  $R(T).  An  example  of  this  is  given  in  [23],  where 

a function  * (T)  is  given  as  exp (1.405-0. 002 7T) . Thus,  the  thermal  and 

K 

irradiation  reduced  times  and  £R  are  now  replaced  respectively  by 
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^ and  ?'  = !*%&&))%(*<*'»&'  (14) 

Utilizing  these  replacements  in  Equations  (10)  and  (13) , and  assuming  that 
the  results  may  be  superposed,  we  obtain 


£u>,a. 


E t.,  o ‘t  r *T 


A U)  A 


*£  I*  6 »j,5 


(15) 


The  material  characterized  by  this  relation  will  be  termed  thermof luxorheologically 
simple.  In  some  cases,  it  may  also  be  necessary  to  generalize  the  thermal 
expansion  and  swelling  terms  through  the  introduction  of  coupled  functions  such 
as  (v,o).  There  appears  to  be  no  experimental  evidence  suggesting  that 

the  thermal  expansion  term  aT  0 be  modified  to  include  a dependence  on  v. 
However,  temperature  may  enter  into  the  swelling  term  as  in  [24],  where  the 
factor  exp  | [ S-Qi j- [ j is  introduced.  In  this  paper,  we  assume  that 


<*T  and 


are  independent  of  0 and  v. 


Converting  Equation  (15)  into  real  time  space,  we  obtain  the  final  result 
£ * - 

c i«i  ° 

(it 

(16) 


-f  + .1,  i Jt‘<  V«  - V*'» 

+|t  j"  y;’(  ?,<«- 


+ o(T  e + o<R  v 

This  result  differs  from  that  given  in  [10]  in  two  respects:  (1)  Two 

physically  different  reduced  time  scales  are  present  rather  than  one,  and 
(2)  The  nonlinear  stress  dependence  is  obtained  by  generalizing  stress  to 
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functions  of  stress  rather  than  by  introducing  stress  into  the  compliance. 
We  now  seek  to  show  that  Equation  (16)  is  consistent  with  a thermodynamic 
development  of  three-dimensional  constitutive  relations. 
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3.  GENERAL  MEMORY  INTEGRAL 

In  the  absence  of  neutron  irradiation,  the  first  law  of  thermodynamics 
for  infinitesimal  theory  (e.g.  see  [15])  is  given  by  the  equation 


-h 


"d  ffij 


‘J  at  aXi 


f ('lr4si r+TirJ 


(17) 


and 


Here,  is  the  strain  tensor,  the  stress  tensor, 

the  divergence  of  heat  efflux  vector  Q^,  p the  mass  density  (assumed 

constant),  G the  Gibbs  free  energy  ( G * ' u +'Ts  * with  °ij 

T as  independent  variables),  s the  entropy  and  u the  Internal  energy. 

When  a metal  is  exposed  to  neutron  flux,  additional  energy  is  transferred 

throughout  the  mass  (by  collision)  primarily  into  two  effects:  (1)  The 

generation  of  local  heating  via  lattice  vibrations,  and  (2)  The  production  of 

point  defects.  The  rate  of  energy  addition  via  the  first  effect  may  be 

expressed  simply  by  the  volume  integral  j^rdv,  where  r is  the  strength  of 

a distributed  internal  heat  source  [25],  In  considering  the  second  effect,  it 

seems  reasonable  to  postulate  that  the  rate  of  energy  added  per  unit  volume 

Cr-f 

in  the  creation  of  point  defects  is  proportional  to  impection  density.  Thus, 
we  may  express  the  rate  of  energy  addition  due  to  the  second  effect  as 


(18) 


where  the  negative  sign  is  taken  for  convenience  of  later  development  and  e 
is  a material  property  (the  conjugate  property  to  n). 

The  global  form  of  the  first  law  of  thermodynamics  is  now  given  as 


• e • • 

P + Q + R - U (19) 

e 

Here,  in  addition  to  the  usual  rate  of  work  (P)  and  time  rate  of  change  of 
total  internal  energy  (U)  terms,  there  is  the  rate  of  heat  influx  term  (Q) 
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as  modified  by  the  addition  of  a heat  source  component, and  the  new  term 
(R)  as  previously  defined  in  Equation  (18).  For  infinitesimal  theory,  the 


P,  Q and  U terms  may  be  written  as 


P * "S 


(20a) 


+ y prdy 

v-  ax;  -V 


(20b) 


u =^f  j ir^+fv 


(20c) 


where  is  the  body  force  and  is  the  velocity.  Combining  Equations 

(18),  (19),  (20)  and  the  definition  of  G in  the  usual  manner,  we  obtain  the 

extented  version  of  local  Equation  (17)  as 

£{J  at  axi  ' at  ' 1 at  5 at  (2i) 

where  now  the  independent  variables  for  G are  o^ , T and  n.  The 
second  law  is  expressed  by  the  usual  Claus lus-Duhem  inequality  with  a heat 
source  term 

SS  9Qi  Qi  aT 

f-r-ti-- m 

Combining  Equations  (21)  and  (22),  and  introducing  v and  0 as 
previously  defined,  we  obtain  for  infinitesimal  theory 


' S at  r at" 


rs  at  t at  e at  20 


This  inequality  forms  the  basis  of  the  development  in  this  section,  where 
now  the  Independent  variables  are  o^  , 8 and  v. 

Proceeding  in  a manner  similar  to  that  in  [15],  as  modified  by  the 
addition  of  the  new  independent  variable  v,  we  express  the  Gibbs  free 


a 
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energy  as 


4, 


f G l <r^fc.t>,0<Xi.*),V(  *.*)]  + G 

where  G*  Is  the  Gibbs  free  energy  function  (per  unit  volume)  due  to  the 

M 

instantaneous  response,  and  G the  Gibbs  free  energy  function  due  to 

M 

memory.  We  postulate  for  G the  sum  of  nonlinear  memory  integrals 

g"- i * 


(24) 


r-( 


(25) 


r<rt 

where  4 m*  , 

-lr) 


are  tensor  functions  of  , 0 and  v,  and 


the  J"  are  time  dependent  material  properties. 

Applying  Leibnitz's  rule  to  Equations  (24)  and  (25),  we  obtain  the  time 
derivative 

3G  . 9C>  90  Xj  VL 

PUT  m Tfc  Ti+  90  dv  a-t 


4 f ^ m 2*  + if£  vL ) ct-t' 

L + IT  3V  9*  J l 9 it-t') 

+ , (-TO  L ii-  JU-mt 


(26) 


We  now  substitute  Equation  (26)  into  inequality  (23)  and  group  terms  to 
obtain 


t ^ + a<j  k 3<rt) )a  '**  at-e  -to  > at 

‘ rS  -99  +£,  o *»  -t')  J 


?e 

rt 


+ ( - c + 2*  *£  2^1 22L 

rt?  ""  '**  O(-fc-t') J-k!=t 

+ £ f‘°  f*  ain  aJlL(^)-77v))  , ^ ae 


(27) 
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Inequality  (27)  must  hold  for  arbitrary  values  of  , d ‘Vjf  and  99/-jt  ; 

this  then  implies  the  constitutive  relations 


Ps 


e » 


•af-t-t'j 


‘I#  ‘I, 

4 +k 


(28a) 


(28b) 

(28c) 


As  a consequence  of  Equation  (28) , inequality  (27)  now  requires 


_A 

where  A 


- JLl  > o 

TV 


is  the  dissipation  function  given  by 


(29) 


<M 


(30) 


Since  Inequality  (29)  must  be  satisfied  for  any  case,  we  may  consider  the 
special  case  of  uniform  temperature  field  =.  o which  then  requires 


A > o 


(31) 


Returning  again  to  arbitrary 


, we  obtain  from  (29)  and  (31) 


fl,  2§_ 

T,  IK 


20 


This  result  is  consistent  with  the  Fourier  heat  conduction  law 


(32) 


GU 


L . 21. 


and  was  obtained  in  a similar  manner  by  Christensen  in  [26]. 


(33) 
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With  the  use  of  Equations  (24),  (25)  and  (28),  coupled  energy  equation 
(21)  may  now  be  rewritten  for  infinitesimal  theory  as 


f r - -§£*  + A -T, rM—  +Z  -^jY*  _ 


(34) 


where  the  term  in  the  brackets  is  the  entropy  s.  A coupled  energy  equation 
for  linearized  infinitesimal  theory  may  be  obtained  by  dropping  the  second 
order  term  A. 

In  the  following  sections,  we  will  show  that  superposition  integral 
(16)  is  in  fact  a special  case  of  Equation  (28  a). 


-17- 


* 


4.  MODIFIED  SUPERPOSITION  INTEGRAL 

In  this  section,  we  obtain  a three-dimensional  constitutive  equation 
containing  superposition  integrals  similar  to  those  in  Equation  (16).  To 
this  end,  we  will  set  in  Equations (28) 


P - M ■ 

f N 

(35a) 

*Tmo 

(35b) 

Cr) 

^.(r) 

; <F-V'i> 

r*  i,  2,  m 

* 

1 

<r) 

• (Fs)tl(  ffij) 

r = M + l,  ...  f SA  ♦ M 

(35c) 

(rt 

r - J,  2,  ••• , M 

1 

r » Mti , ••*,  M ♦ N 

(35d) 

It 

H 

Or 

G + £*  oij  + 

fse  + ev  -»-^CjIu<ra<r«*  +5-^0' 

+ j * (dr\j  CTy  9 + 


+ Zljj^0)f\Fy)l”itr!j)dor~  + Z \fJ2(o$)  day 


(35e) 


Here  (F_)  ^ and  (F_)  (r) 

1 0,11  R 1,1111  are  functions  of  stress  only,  is  the 

linear  elastic  compliance,  (a.j.)^  the  thermal  expansion  coefficient,  Cq 
the  specific  heat  at  constant  stress,  G,  , s and  e are  initial  values, 
and  p,  h and  (Og)^  are  new  material  constants. 

Employing  Equations  (35a- e)  in  Equations  (28)  we  obtain 


CjM  or*  * 9 + 5 A j^(  Vi)  - at' 

N t IT)  Jrf.l 

+ z I,  < J. <*  > ' ?« <*'> ) at ' 


(36a) 
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PS  - f>S  + ( 4r)ii  <Ty  + + 


hx? 


(36b) 


e = e * <«>(*)-  (Ty  + he  + 


(36c) 


where  we  have  used  an  integration  by  parts  and  the  conditions  (0~)« 
£^(0  )m0.  Inserting  the  entropy  (36b)  Into  energy  equation  (34),  we  also 
obtain 


r_-Z£itA.fCJ».-Wi2g.-Trhl2-.o 


(37) 


Equation  (36a)  is  clearly  the  desired  extension  of  constitutlve„equatlon 
(16)  to  three  dimensions.  Entropy  expression  (36b)  is  the  same  as  in  linear 
thermoelasticity,  with  the  addition  of  the  structure  dependent  term  hv. 
Finally,  if  we  drop  the  second  order  dissipation  term  from  coupled  energy 
equation  (37),  we  obtain  the  result  of  linear  thermoelasticity  supplemented 

by  the  two  new  energy  terms  and  -Trh|^-' 

For  an  isotropic  material,  Fourier  heat  conduction  equation  (33) 

reduces  to 


Qi 


and  then  energy  equation  (37)  becomes 


(38) 


fr  + k-^f  + A- 


(39) 


Finally,  the  stress-strain  relation  (36a)  is  written  for  an  Isotropic 
material  as 
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+ Z 

*»i 


tvt  J . , /ir/ 

f»i  • 

’Z  iHVv1'11 * 5 l<1i(Wu-«t,KI«fi?^*'^  !^tilLi  M'  («» 

r*»  • 


Here,  vg  and  E are  the  elastic  Poisson's  ratio  and  Young's  modulus 
respectively,  and  the  subscripts  I and  II  designate  shear  and  bulk  creep 
compliance  functions  respectively.  In  the  next  section  we  discussed  the 
power  law  form  of  Equation  (40). 
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5.  ISOTROPIC  STRESS  POWER  LAW 


In  [27],  a nonlinear  constitutive  equation  of  the  isotropic  power 
law  type  was  developed  for  a three-dimensional  viscoelastic  material  with 
instantaneous  and  time  dependent  compressibility.  The  development  was 
based  on  the  assumption  that  creep  potentials  may  be  obtained  by  raising 
the  elastic  quadratic  form  to  integer  powers,  and  furthermore  the  creep 
compliances  used  were  for  the  generalized  Kelvin  model.  In  [15],  Chang 
et  al.  equated  the  Integrated  function  in  the  Gibbs  free  energy  (in  our 
case  gjj"(F^)^  or  (Fr)j^  in  Equations  (25,  36a,  40))  to  the  derivative 
(■^ — ) of  the  above  mentioned  potentials,  and  obtained  the  same  single 
integral  representation.  In  this  section  we  shall  also  use  these  same  forms 
for  the  functions  (F^).^  and  (FR) in  Equation  (40),  and  furthermore  we 
will  also  utilize  generalized  Kelvin  compliances. 

Proceeding  as  described  above,  we  obtain  the  isotropic  stress  power 

law 


V 


■v» 


>*  - Vi M 

* s^i,*jVC  r,  - 7&-I,  VK 


(41) 
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In  the  above,  is  the  first  invariant  of  stress  and  J0«l/2s . .s, . 

1 ‘ ij  ij 

is  the  second  invariant  of  the  stress  deviator  s . *0 . .-l/3a, . 6,. . . Sub- 

ij  ij  kk  ij 

scripts  Ts,  Tt,  Rs  and  Rt  designate  thermal  steady,  thermal  transient, 
irradiation  steady  and  irradiation  transient  creep  quantities  respectively. 
Accordingly,  the  various  creep  compliances  are  defined  as 

JfS  3 


7*(<>  ^ 3 


> I - e*Fl-CV'>) 


k*  2,  vl 

k - 2,*”,  N 


(42) 


in  conformity  with  Equations  (1)  and  (11).  The  powers  m^,  m^^ , m^  and 


.00 


mRt  are  all  restricted  to  integers  0,1,2,...;  the  constants  v , v**' , \>„k) 

i 8 It  Kt  9 

are  Poisson's  coefficients  for  the  various  creep  components;  and  the  C,pg, 

CTt  ’ CRs  an<*  CRt^  *re  creeP  constants  which  multiply  the  various  compliances. 

The  time  dependent  compressibility  is  displayed  by  forming  the 
contraction  of  Equation  (41) 

tu  - -TT*-  + 

A)  *' 

tv  v»  ■ *' 


rKt  • 

+ 3 o(„  + 3 o lr6 


(43) 
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This  equation  is  useful  in  analyzing  pressure  test  data. 

For  the  one-dimensional  creep  test,  a ma  H(t)  where  a is  constant 

xx  o o 

and  H(t)  is  the  unit  step  function,  and  the  flux  and  temperature  are  also 
kept  constant  at  and  Tq.  Equation  (41)  then  yields 

_ <k)  w 


N o 

+w  7s ( l - ]}  Htt)  + 0^5,+  ol„  4 

Jk) 

= £«  = * [ >2  y + +2 ; (^t)  [ I - exp(- 1 aJ^t  )] 

4V«*(^  W *£$^1 1 - e*p(-ti%%>)}  Hrt)  +M+ 


(44a) 


(44b) 


In  Equations (44)  we  have  eliminated  reduced  times  nt  and  by  means  of 
Equation  (14)  with  constant  and  and  we  have  also  Introduced 


nT=2mH+l 


J»)  ik) 


^Va-yo  ( 1 1 I ) 

y(z 

* r *» 


' r (/  (*)  . 
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The  role  of  the  various  Poisson's  coefficients  becomes  clear  when  we 

compare  Equation  (44b)  with  (44a).  In  [28J  it  is  suggested  that  v_  is  close 

Ts 

(k) 

to  1/2  while  v^t  lie  somewhere  between  vg  and  1/2.  Although  some  Information 
on  v is  given  in  [19],  there  appears  to  be  no  data  available  on  v„  ^ . On 
comparing  Equation  (44a)  with  Equations  (1)  and  (11),  we  see  that  we  have 
essentially  obtained  our  initial  equations.  As  remarked  earlier,  irradiation 
creep  is  usually  taken  as  linear  in  stress.  Thus  we  may,  with  little  loss 
of  accuracy,  set  nR«qR  »1  in  Equations  (44).  Finally,  one  may  readily 
specialize  energy  equation  (39)  for  the  Isotropic  power  law  by  evaluating 
dissipation  function  (30)  for  this  material,  but  as  previously  noted  it  is 
common  practice  to  simply  drop  the  dissipation  function  from  the  energy  equation. 
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6.  SUMMARY 


One-dimensional  constitutive  equations  valid  for  thermal  and  irradiation 
creep  at  constant  stress  were  generalized  to  variable  stress  by  utilizing 
the  strain  hardening  technique  to  obtain  modified  superposition  integrals. 
Thermal  creep  terms  valid  for  variable  temperature  were  obtained  first  by 
assuming  a thermorheologically  simple  material, accordingly  introducing 
a thermal  reduced  time  to  account  for  temperature  dependence  of  the  material 
parameters.  Similarly,  irradiation  creep Avalid  for  variable  flux  were 
obtained  by  assuming  a fluxorheologically  simple  material,  and  introducing  an 
irradiation  reduced  time  dependent  on  the  A rate  of  change  of  the  imperfection 
density.  An  expression  for  thermoirradiation  creep  (i.e.,  temperature 
and  flux  both  high)  was  then  obtained  by  superposition  of  the  previous  thermal 
and  irradiation  creep  terms,  as  modified  by  a redefinition  of  the  two 
reduced  time  scales  to  account  for  coupling.  Linear  elastic,  thermal 
expansion  and  swelling  terms  were  also  included  in  the  one-dimensional 
constitutive  equation. 

A systematic  development  based  on  the  first  and  second  laws  of 
thermodynamics  then  enabled  us  to  obtain  general  memory  Integrals  which 
contain  as  a special  case  those  previously  obtained  in  the  paper  by 
physical  argument.  The  conventional  forms  of  the  infinitesimal  lavs  of 
thermodynamics  were  first  modified  to  account  for  the  effect  of  neutron 
flux  on  local  heating  and  point  defect  production.  A Gibbs  free  energy 
was  assumed  as  the  sum  of  an  Instantaneous  function  component  and  a memory 
functional  component  with  Btress,  temperature  increment  and  defect  density 
increment  as  independent  variables.  As  a result,  three-dimensional  an- 
isotropic constitutive  equations  were  obtained  for  the  strain,  entropy  and 
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a material  property  conjugate  to  the  defect  density;  the  coupled  energy 
equation,  dissipation  function  and  heat  conduction  law  were  also  obtained. 

Modified  superposition  integrals  similar  to  those  previously  developed 
were  then  obtained  by  substituting  special  functions  into  the  results  of 
the  thermodynamic  analysis.  Finally,  an  Isotropic  stress  power  law  was 
obtained  through  the  use  of  quadratic  forms  raised  to  integer  powers,  and 
was  then  examined  for  three-dimensional  time  dependent  compressibility  and 
one-dimensional  Poisson  contraction. 

The  validity  of  the  analogy  between  thermally  induced  and  irradiation 
induced  creep  as  presented  here  must  of  course  be  tested  against  further 
experimental  evidence.  However,  we  have  conclusively  shown  that  the  analogy 
as  deduced  from  simple  one-dimensional  data  is  consistent  with  a more 
general  thermodynamic  development.  We  also  point  out  that  the  analogy 
should  prove  useful  in  reducing  the  complexity  of  future  testing  programs. 
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* Employing  an  analogy  between  thermally  induced  and  irradiation  induced 
creep,  physical  arguments  are  used  first  to  deduce  a one-dimensional 
constitutive  relation  for  metals  under  stress  in  a high  temperature  and  high 
neutron  flux  field.  This  constitutive  relation  contains  modified  superposition 
Integrals  in  which  the  temperature  and  flux  dependence  of  the  material  parameters 
is  included  via  the  use  of  two  reduced  time  scales;  linear  elastic,  thermal 
expansion  and  swelling  terms  are  also  included.  A systematic  development  based 
on  thermodynamics,  with  the  stress,  temperature  increment  and  defect  density 
increment  as  independent  variables  in  the  Bibbs  free  energy,  is  then  employed 
to  obtain  general  three-dimensional  memory  integrals  for  strain;  the  entropy  and 
coupled  energy  equation  are  also  obtained.  Modified  superposition  integrals 
similar  to  these  previously  obtained  by  physical  argument  are  then  obtained  by 
substituting  special  functions  into  the  results  of  the  thermodynamic  analysis, 
and  the  special  case  of  an  isotropic  stress  power  law  is  examined  in  detail. 


DD  ,'",“..1473  " 


/N  AIM 


it  AY  I III  I 


UNCLASSIFIED 


I'AKV  i - GOVERNMENT 


Anni ni  strrt.ivc  A Liaison  Activities 

Thief  of  Naval  Research 
Jjcj-artr.rnt  of  the  Navy 
Arlington,  Virginia  22217 
Attn:  Code  1*39  (2) 

li71 

Di  rector 

ONR  Branch  Office 
1*95  Sumner  Street 
Boston,  Massachusetts  02210 

Director 

OUt  Brunch  Office 

219  G.  Dearborn  Street 

Chicago,  Illinois  6060U 

Director 

Koval  Research  Laboratory 
Attn:  Library,  Code  2029  (ONRL) 
Washington,  D.C.  20390  (6) 

Co,:imandi ng  Officer 
ON  It  branch  Office 
207  best  2l|th  Street 
Now  York,  New  York  10011 

Director 

ONR  Branch  Office 
1030  E.  Green  Street 
Pasadena,  California  91101 

U.S.  Naval  Research  Laboratory 
Attn:  Technical  Information  Div. 

Washington,  D.C.  20390  (6) 

Defense  Documentation  Center 
Cameron  Station 

Alexandria,  Virginia  22311*  (20) 

Army 

Commanding  Officer 

U.S.  Army  Research  Office  Durham 

Attn:  Kr.  J.J.  Murray 

CRD-AA-IP 

Box  GK,  Puke  Station 
Durham,  North  Carolina  27706 


Command  i ng  Offi  cer 
AMXKR-ATL 

Attn:  Kr.  J.  Bluhm 

U.S.  Army  Materials  Res.  Agcy. 

Watertown,  Massachusetts  02172 

"Watervlict  Arsenal 
MAGGS  Research  Center 
Watervlict,  New  York  12189 
Attn:  Director  of  Research 

Restonc  Scientific  TnJo.  Center 
Chief,  Document  Section 
U.S.  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35809 

Arm y R & D Center 

Fort  Bel voir,  Virginia  22060 

Technical  Library 
Aberdeen  Proving  Ground 
Aberdeen,  Maryland  21005 


Commanding  Officer  and  Director 
Naval  Ship  Research  & Development.  Cer. 
Washington,  D.C.  20007 
Attn:  Code  0l*2  (Tech.  Lib.  Br.) 

700  (Struc.  Kech.  Lab.) 
720 

725 

727 

012.2  (Dr.  W.J.  Sctte) 

Naval  Weapons  Laboratory 
Dahlgrcn,  Virginia  221,1*8 

Naval  Research  Laboratory 
Washington,  D.C.  20390 
Attn:  Code  81:00 

81*10 
81*30 
81,1*0 
6300 
. 6305 
638O 

Undersea  Explosion  Research  Div: 

Naval  Ship  R Sr  D Center 
Norfolk  Naval  Shipyard 
Portsmouth,  Virginia  23709 
Attn:  Dr.  Schauer 

Code  7C0 


Naval  Ship  K ft  D Center 

Annnpol in  Divi  u Lon 

Annapolis,  Maryland  21U02 

Attn:  Code  AOOO,  Mr.  W.L.  ’..Viliams 

Technical  Library 

Naval  Underwater  Weapons  C'  er 

Pasadena  Annex 

3202  E.  Foothill  Blvd. 

Pasadena,  California  9110 < 

U.S.  Naval  Weapons  Center 
China  Lake,  California  9 3997 
Attn:  Code  1(920  Mr.  Ken  Bisehel 

U.S.  Naval  Ordnance  Laboratory 
Mechanics  Division 
RrT)  1,  White  Oak 
Silver  Spring,  Maryland  20910 

U.S.  Naval  Ordnance  Laboratory 
Attn:  Mr.  H.A.  Perry,  Jr. 

Non -Motalii c Materials  Division 
Silver  Spring,  Maryland  20910 

Technical  Director 

U.S.  Naval  Undersea  K ft  D Center 

San  Diego,  California  92132 

Supervisor  of  Shipbuilding 
U.S.  Navy 

Newport  News,  Virginia  2360 7 

Technical  Director 

Mare  Island  Naval  Shipyard 

Vallejo,  California  9U992 

U.S.  Naval  Ordnance  Station 
Attn:  Mr.  Garct  Borns Lvi*, 

Research  ft  Development  jjj.v 
Indian  Head,  Maryland  ?06li0 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington,  D.C.  20390 
Attn:  Code  0P-07T 

Deep  Submergence  Systems 
Naval  Ship  Systems  Command 
Code  39922 

Department  of  the  Navy 
Washington,  D.C.  PC;"' 

Attn:  Chief  Scienta." 


Director,  Aero  Mechanics 
Naval  Air  Development  Center 
Johnsvillc 

Warminster,  Pennsylvania  1897U 

Naval.  Air  Systems  Command 
Dept,  of  the  Navy 
Washington,  D.C.  20360 
Attn:  NAIR  320  Aero,  ft  Structures 

9320  Structures 
60Jt  Tech.  Library 
92031F  Materials 

Naval  Facilities  Engineering  ( 

Command 

Department  of  the  Navy 

Washington,  D.C.  20360  , 

Attn:  NE'AC  03  Res.  ft  Development 

OU  Engineering  ft  Design 
lhllij  Tech.  Library 

Naval  Ship  Systems  Command 
Dept,  of  the  Navy 
Washington,  D.C.  20360 
Attn:  NSHTP  031  Ch.  Scientists  for  P 

03U2  Ship  Mats,  ft  Struct j 
209'2  Tech.  Library 

Naval  Ship  Engineering  Center 
Prince  George  Plaza 
Hyattsvillc,  Maryland  20782 
Attn:  NSEC  6100  'hip  Sys.  F.ngr.  ft  D< 

6102C  Computer-Aided  Shi 

6109  Ship  Protection 

6110  Sfiio  Concept  Design 
6120  Hull  Div. 

612UD  Hull  Div. 

6128  Surface  Ship  Struct 

6129  Submarine  Struct. 

Naval  Ordnance  Systems  Command 
l)cpt.  of  the  Navy 
Washington,  D.C.  20360 
Attn:  NORD  03  Res.  ft  Technology 

039  Weapons  Dynamics 
9132  Tech.  Library 

Engineering  Department 
U.S.  Naval  Academy 
Annapolis,  Maryland  21U02 


-4- 


Air  • 

Conniandcr  WADD 

Wrignt-Patt«  rson  Air  Force  Base 
Dayton,  Ohio  h51;33 
Attn : Code  WWKMDD 

AFFDL  (FDDD) 

Structures  Division 
AFLC  (MCSEA) 

Code  YMRC 
AFIIL  (maam) 

Ccrannder 

Chief,  Applied  Mechanics  Group 
U.S.  Air  Force  Inst,  of  Tech. 
Wright-Palterson  Air  Force  Base 
Dayton,  Ohio  hlJh33 

Chief,  Civil  Engineering  Branch 
VLRC,  Research  Division 
Air  Force  Weapons  Laboratory 
Kirfland  AFB,  Kev;  Mexico  87117 

Air  Force  Office,  of  Scientific  Res. 

U;00  Wilson  Diva. 

Arlington,  Virginia  22209 
Attn:  Mechs.  Div. 

t.sf.A 

Structures  Research  Divi  sion 
Rational  Aeronautics  Ac  Space  Admin. 

.]  Langley  Research  Center 
i hartley  Station 

• Bamptori,  Virginia  233^5 

• Attn:  Mr.  R.R.  Hcldenf els.  Chief 

National  Aeronautic  & Space  Admin. 
Associate  Anninistrator  for  Advanced 
Research  ft  Technology 
Washington,  D.C.  ?0fJI\6 

Scientific  & Tech.  Info.  Facility 
RASA  Representative  (S-AK/iJi.) 

! P.O.  Box  5700 

Bcthosda,  Maryland  2001^ 

Other  Covermont  Activities 

• Technical  Di  rector 

Marine  Corpr.  Development  At  Educ.  Command 
Ounntico,  Virginia  221  'jh 


Director 

National  Bureau  of  Standards 
Washington,  D.C.  2023U 
Attn:  Mr.  3.L.  Wilson,  EM  219 

National  Science  Foundation 
Engineering  Division 
Washington,  D.C.  20550 

Director 

STDS 

Defense  Atomic  Support  Agency 
Washington,  D.C.  20350 

Commander  Field  Command 
Defense  Atomic  Support  Agency 
Sandia  Base 

Albuquerque,  New  Mexico  87115 

Chief,  Defense  Atomic  Support  Agcy. 
Blast  At  Shock  Division 
The  Pentagon 
Washington,  D.C.  20301 

Director  Defense  Research  A:  Engr. 

Techni cal  Library 

Room  30-128 

The  Pentagon 

Washington,  D.C.  20301 

Chief,  Airfra.de  A:  Equipment 
FS-120 

Office  of  Flight  Standards 
Federal  Aviation  Agency 
Washington,  D.C.  20553 

Chief  of  Research  and  Development 
Maritime  Administration 
Washington,  D.C.  20235 


Mr.  Milton  Shaw,  Director 
Div.  of  Reactor  Dcvolo.  & Technology 
Atonic  Fiiergy  Commission 
Gomantovm,  Maryland  20?67 

Ship  Hull  Research  Committee 
National  Research  Council 
National  Academy  of  Sciences 
2101  Constitution  Avenue 
Washington,  D.C.  20U18 
Attn:  Mr.  A.R.  Lytle 


4 


PAirr  2 - or  infractors  and  other 

TECHNICAL  COLLABORATORS 
Dll  I Vcr.Tj  tif;S 

Professor  J.R.  Rice 
Division  of  Engineering 
Brovin  University 
Providence,  Rhode  Island  02912 

Dr.  J.  Tinsley  Oden 
Dept,  of  Engr.  Kochs. 

University  of  Alabama 
Huntsville,  Alabama  35>0Oli 

Professor  K.S.  Rivlin 

Center  for  the  Application  of  Mathematics 

Lehigh  University 

Bethlehem,  Pennsylvania  18015 

Professor  Julius  1-Iiklowitz 
Division  of  Kngr.  ft  Applied  Sciences 
California  Institute  of  Tec'nnology 
Pasadena,  California  91109 

Professor  George  Sih 
Department  of  Mechanics 
Z .chi  j;h  Uni. v *•  r sity 
Beti-il  eh  cm,  Pennsylvania  18015 

Dr.  Harold  Liebowitz , Doan 
School  of  Engrg.  & Applied  Science 
George  Washington  University 
725  23rd  Suroet 
’Washington,  D.  C.  20006 

iVofcssor  Eli  Sternberg 
Division  of  Engr.  L Applied  Sciences 
California  Institute  of  Technology 
Pasadena,  California  91109 

Professor  Burt  Paul 

University  of  Pennsylvania 

Tovmc  School  of  Civil  and  Mechanical  lingrg. 

Room  113  Towne  Building 


?20 


33rd  Street 


Philadelphia,  Penna.  1910[i 

Professor  S.B.  Dong 
University  of  California 
Department  of  Mechanics 
Los  Angol.cs,  California  90021* 


Professor  Paul  M.  K'aghdi 
Div.  of  Applied  Mechanics 
Etchcvcrry  Hall 
University  of  California 
Berkeley,  California  91*720 

Professor  W.  Nachbar 
University  of  California 
Department  of  Aerospace  L Koch,  Engrg 
La  Jolla,  California  92037 

Professor  J.  Baltrukonis 
Mechanics  Division 
The  Catholic  Univ.  of  America 
Washington,  D.C.  20017 

Professor  A.J.  IXirclli 
Mechanics  Divi sion 
The  Catholic  Un.iv.  of  America 
Washington,  D.C.  20017 

Professor  H.H.  Bleich 
Department  of Civil  Engrg. 

Columbia  University 
Amsterdam  ft  120th  Street 
New  York,  Keu  York  10027 


Professor  R.D.  Mindlin 
Department  of  Civil  Engrg. 

Columbia  University 

S.W.  Kudd  Building 

New  York,  New  York  1002? 

Professor  A.M.  Ereudenthal . 

Oeorge  Washington  University 
School  of  Fhgrg.  & Applied  Science 
Washington,  D.C.  20006 

Professor  B.A.  Boley 
Department  of  Thoor.  & Appl.  Mcch. 
Cornell  University 
Ithaca,  New  York  11*850 

Professor  P.Q.  Hodge 
Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago,  Illinois  60616 

Dr.  D.C.  Drucker 
Doan  of  hhgineoring 
University  of  Illinois 
Urbana,  Illinois  61801 


-it- 


•fcTTvwt  * Trr^-^r*!  7.pv» 


Professor  f.'.K.  N<  wr.-nr/ 

Oc-pi- . ofCivil  Engineering 
University  of  Illinois 
Urbana,  Illinois  61801 

Professor  Jnr.es  Mar 
Massachusetts  Inst,  of  Tech. 

Rn.  33-318 

Dept,  of  Aerospace  U Astro. 

77  Massachusetts  Avc. 

Cambridge,  Kassachusettc  02139 

Library  (Code  0j”J>) 

li.S.  Naval  Postgraduate  School 

Monterey,  California  939UO 

Dr.  Francis  Cor.zarclli 
Div.  of  Interdisciplinary 
Studies  and  P.escarch 
School  of  fckiginocring 
State  Univ.  of  N.Y. 

Buffalo,  New  York  lii2lU 

Professor  K. A.  Douglas 
Dept,  of  Engrg.  Kochs. 

North  Carolina  State  University 
Raleigh,  North  Carolina  27607 

Dr.  George  Herrmann 
Stanford  University 
Department  of  Applied  Mechanics 
Stanford,  California  9lj 30^’ 

Professor  J.D.  Achonbach 
Technological  Institute 
Northwestern  University 
Evanston,  Illinois  60201 

Director,  Ordnance  Research  Lab. 
Pennsylvania  State  University 
P.0.  Box  30 

State  College,  Pennsylvania  16801 

Professor  J.  Kcsapner 

Dept,  of  Aero.  Engrg.  & Applied  Mech. 

Polytechnic  Institute  of  Brooklyn 

333  Jay  Street 

Brooklyn,  New  York  11201 

Professor  J.  Klosncr 
Polytechnic  Institute  of  Brooklyn 
333  Jay  Street 
Brooklyn,  New  York  13.201 


Professor  A.C.  Eringen 

Dept,  of  Aerospace  & Kcch.  Sciences 

Princeton  University 

Princeton,  New  Jersey  085UO 

Dr.  S.L.  Koh 

School  of  Aero.,  Astro.  & Engr.  Sc. 
Purdue  University 
I.afayett,  Indiana  1(7907 

Professor  R.A.  Schapcry 
Civil  Jhg:i  neering  Department 
Texas  A & K University 
College  Station,  Texas  7781(0 

Professor  K.H.  Lee 
Div.  of  Engrg.  Mechanics 
Stanford  University 
Stanford,  California  9l(30S> 

Dr.  Nicholas  J..Hoff 
Dept,  of  Aero,  & Astro. 

Stanford  University 
Stanford,  California  9h30$ 

Professor  Max  Anliker 
Dept.  of  Aero.  Sc  Astro. 

Stanford  University 
Stanford,  California  9U303” 

Professor  Chi. -Chang  Chao 
Div.  of  Erigr.  Mechanics 
Stanford  University 
Stanford,  California  9U3 0J> 

Professor  H.W.  Liu 

Dept,  of  Chemical  Engr.  & Metal. 

Syracuse  University 

Syracuse,  New  York  13210 

Professor  S.  Bodnor 
Tcchnion  R & D Foundation 
Haifa,  Israel 

Dr.  S.  Dhawan,  Director 
Indian  Institute  of  Science 
Bangalore,  India 

Professor  Tsuyoshi  Hayashi 
Department  of  Aeronautics 
Faculty  of  Engineering 
University  of  Tokyo 
BUNK YD -KU 
Tokyo,  Japan 


Professor  J.F..  Fil  gr 
Dept.  of  Civil  Knj  . <>•■ 
University  of  Utah 
Salt  Lake  City,  U«  h 


Professor  R.J.H.  P ‘‘  l 
Chairman,  Aeronaut  ca.  gr 
20V  Guggenheim  Hall 
University  of  Washington 
Scattic,  Washington  93105 


Dept. 


Professor  Albert  S.  Kobayusni 
Dept,  of  Mechanical  Engr. 
University  of  Washington 
Seattle,  Washington  98105 


Professor  G.R.  Irwin 
Dept . of  Mech.  Engrg. 
Lehigh  University 
Bethlehem,  Pennsylvania 


l801j; 


Dr.  Daniel  Frederick 
Dept,  of  Engr.  Keens. 
Virginia  Polytechnic  Inst. 
Blacksburg,  Virginia  2U061 


Professor  Lambert  Tall 


Lehigii  University 
Department  of  Civil  Engrg. 
Bethlehem,  Pennsylvania  18015 


Professor  K.P.  Wnuk 
South  Dakota  State  University 
Department  of  Mechanical  Engineering 
Brookings,  South  Dakota  57006 


Professor  Norman  Jcr.es 
Massachusetts  Insitute  of  Technology 
Dept,  of  Naval  Architecture  & Marine  Ehgrg . 
Cambridge,  Massachusetts  02139 


Professor  Pedro  V.  Karcal 
Brovrn  University 
Division  of  Engineering 
Providence,  Rhode  Island  02912 


Professor  Werner  Goldsmith 
Department  of  Mechanical.  Engineering 
Division  of  Applied  Mechanics 
University  of  California 
Berkeley,  Cal.ifornia  91:720 


Professor  R.B.  Testa 
Dept,  of  Civil  Engineer! ng 
Colombia  University 
S.W.  Xudd  Bldg. 

New  York,  Mew  York  10027 


-6- 


Dr.  Y.  Weitsnan 
Dept,  of  ICngrg . Sciences 
Tcl-Aviv  University 
R. -mat -Aviv 
Tcl-Aviv,  Israel 

Professor  V'.D.  Pilkey 
Dept,  of  Aerospace  Kngrg. 
University  of  Virginia 
Charlottesville,  Virginia  22903 

Professor  W.  Prnger 
Division  of  Engineering 
Broym  University 
Providence,  Rhode  Island  02912 

Industry  andRoscarch  Institutes 

Mr.  Carl  E.  Hartbowcr 
Dept.  1:620,  Bldg.  2019  A2 
A ero  j ct -General  C orporation 
P.0.  Box  19U7 

Sacramento,  California  95809 

Library  Services  Department 
Report  Section,  Bldg,  lli-l): 

Argonne  National  Laboratory 
9h00  S.  Cass  Avenue 
Argonne,  Illinois  60lth0 

Dr.  F.R.  Schwarzl 
Central  Laboratory  T.N.0. 

Schoennak  erstraat  97 
Delft,  The  Netherlands 

Dr.  Wendt 

Valley  Forge  Space  Technology  Cen. 
General  Electric  Co. 

Valley  Forge,  Pennslyvania  10):6l 

Library  Newport  News  Shipbuilding 
tf  Dry  Dock  Company 
Newport  News,  Virginia  23607 

Director 

Ship  Research  Institute 
Ministry  of  Transportation 
700,  SHINKAWA 
Mitaka 
Tokyo,  Japan 

Dr.  II. N.  Abramson 
Southwest  Research  Institute 
8500  Culebra  Road 
San  Antonio,  Texas  78206 


•N 


Dr.  R.C.  DeHart 
Southwest  Research  Institute 
65/CX)  Culcbra  Road 
San  Antonio,  Texas  78206 

Mr.  Roger  Weiss 

High  T«np.  Structurs.  & Materials 
Applied  Physics  Lab. 

6623  Georgia  Ave. 

Silver  Spring,  Maryland  20910 

Mr.  K.C.  Francis,  Head 
Mech.  Props.  Kval. 

United  Technology  Center  ' 

Sunnyvale,  California  94088  , 

Mr.  C.N.  Robinson 
Atlantic  Research  Corp. 

Shirley  Highway  at  Edsall  Road 
Alexandria,  Virginia  222314 

Mr.  P.C.  Dump 
Aeromechanics  Dept.,  74-43 
LocJd  iced-Calif  ornia  Co. 

Burbank,  California  91603 

Mr.  P.  VJilson 
Litton  Systems,  Inc. 

Atf?P,  Dept.  400 
El  Segundo 

9920  V.  Jefferson  Blvd. 

Culver  City,  California  90230 

T)r.  Kevin  J.  Forsberg,  Head 
Solid  Mechanics 
Orgn  5)2-20,  Fldg.  20$ 

Lockheed  Palo  Alto  Research  Lab. 

Palo  Alto,  California  94302 

Dr . K.M.Q.  Roren 
Head,  Research  Department 
Dot  Norske  Veritas 
Post  Pox  6060 
Oslo,  Norway 

Dr.  Andrew  F.  Conn 
Hydronautics,  Incorporated 
Pindcll  School  Road,  Howard  County 
Laurel,  Maryland  20810 


-7- 


